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Abstract 

In this paper, we theoretically prove a set of fundamental conditions pertaining discrete velocity sets and corre- 
sponding weights. These conditions provide sufficient conditions for a priori formulation of lattice Boltzmann models 
that automatically admit correct hydrodynamic moments up to any given A''-th order. 
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1. Introduction 

Lattice Boltzmann methods (LBM) has been 
recognized as an advantageous numerical method 
for performing efficient computational fluid dynam- 
ics |1I2) . Not only it offers a new way of describing 
macroscopic fluid physics, but also it has become a 
practical computational tool and already has been 
making substantial impact in real world engineering 
applications ^3 . Furthermore, according to a more 
recent interpretation, LBM models are special dis- 
crete approximations to the continuum Boltzmann 
kinetic equation |4I5| . Owing to such an underlying 
kinetic theory origin, LBM is expected to contain a 
wider range of fluid flow physics than the conven- 
tional hydrodynamic fluid descriptions [6 7 8 10. 
The latter, such as the Euler or the Navier-Stokes 
equation, rely on various "theoretical" closure ap- 
proximations for the non-equilibrium effects that 
are problematic when deviations from local thermo- 
dynamic equilibrium are no longer considered small. 
In addition, due to the fact that the fundamental 
turbulence modeling is built upon an analogy to 
regular fluid flows at finite Knudsen numbers, a 
kinetic theory representation is argued to be more 
suitable than the classical modeling approach of 



modified Navier-Stokes equations [55]. However, 
how much the original range in kinetic theory can 
be retained depends on the order of accuracy in 
the LBM models used. Indeed, it has been shown 
that certain key physical effects beyond the Navier- 
Stokes equations can be accurately captured using 
higher order LBM models [5l9j . 

There have been extensive studies in LBM for 
more than a decade. However, popularly known 
LBM models are only accurate in the Navier-Stokes 
hydrodynamic regime (c.f., [12ll4j ). That is, physics 
higher than the Navier-Stokes order is contami- 
nated by numerical artifacts in these LBM models. 
Furthermore, there has not been progress in system- 
atically deriving higher order accurate LBM models 
until recently [5]. Originated from the framework 
of the so called Lattice Gas Automata |I5III| . the 
conventional approach to formulating LBM mod- 
els is based on a so called "top down" procedure. 
That is, giving a macroscopic equation such as 
the Euler or the Navier-Stokes equation, an LBM 
model may be constructed via an inverse Chapman- 
Enskog process and a post-priori parameter match- 
ing along with various subsequent "corrections" 
(cf., |13|14|16|17j ). But more fundamentally, be- 
cause such an approach relies on the availability of 
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macroscopic descriptions, it encounters an intrinsic 
difRculty in extending physics beyond the original 
macroscopic equations. It is weh known that there is 
no well established and reliable macroscopic equa- 
tion for deeper non-equilibrium physics beyond the 
Navier-Stokes regime. 

One can theoretically show that the level of non- 
equilibrium physics is directly associated to the 
hydrodynamic moments '5'^. Specifically, from the 
representation of the Chapman-Enskog expansion, 
there exists an apparent hierarchical relationship 
among hydrodynamic moments at various non- 
equilibrium levels. That is, n-th order hydrody- 
namic moments at m-th non-equilibrium level are 
related to the (n-l-l)-th order moments at (m— l)-th 
non-equilibrium level. Carrying out this hierarchy 
all the way, we see that in order to ensure the n-th 
moment physics at the m-th non-equilibrium level, 
it requires the equilibrium moments of (m -|- n)-th 
order to be accurate. In other words, the higher or- 
der equilibrium hydrodynamic moments captured 
accurately, the wider range of non-equilibrium 
physics can be described. Indeed, the popularly 
known lattice Boltzmann models are only accu- 
rate up to the second order equilibrium moment 
(i.e., the equilibrium momentum flux tensor). As 
a result, these models only give an approximately 
correct "level- 1 non-equilibrium momentum flux. 
This is a reason why the conventional LBM models 
are only applicable to the Navier-Stokes (Newto- 
nian) fluid physics in low Mach number isothermal 
situations |11|12|13|14] . 

Based on the above, we see that the essential re- 
quirement for accurately capturing a wider range of 
physics is directly related to achieving equilibrium 
hydrodynamic moments to higher orders. Once the 
higher order moments are accurately realized, the 
resulting hydrodynamic equations such as the Eu- 
ler, the Navier-Stokes and beyond are automatically 
attained. This is accomplished without the conven- 
tional post-priori procedure. As shown in this paper, 
the above requirement dictates a set of fundamental 
conditions on the supporting lattice velocity basis in 
LBM. That is, given an A^-th order moment accu- 
racy requirement, the set of fundamental conditions 
automatically defines the choice of a discrete lattice 
velocity set and its corresponding weights for such 
a purpose. 

In this paper, we theoretically derive this set of 
fundamental conditions for LBM models of A^-th or- 
der. We prove how the correct hydrodynamic mo- 
ments up to the corresponding order are realized 



once the conditions are satisfied. 

2. Achieving Correct Hydrodynamic 
Moments via Discrete Velocities 

According to the standard continuum Boltzmann 
kinetic theory, an n-th order equilibrium hydrody- 
namic moment tensor in Z3-dimension is defined as 



Af(")(a;,t) = J d^ccc^r'^{x,c, t) 



(1) 



Equivalently, it can be expressed in a Cartesian com- 
ponent form as follows. 



^c,J-'^{x,c,t) (2) 



where subscripts ii, Z2, . . . , Zn are Cartesian compo- 
nent indices. Ci is the i-th Cartesian component of 
the microscopic particle velocity c. The equilibrium 
distribution has the standard Maxwell-Boltzmann 
form. 



n{x,c,t) = 



p{x,t) 



X exp 



(c — u{x, t)Y 



(3) 



where the macroscopic density, fluid velocity, and 
temperature are defined, respectively 

p{x,t)= j d''cr'^{x,c,t) 

pu{x,t) = j d^ccr'J{x,c,t) 
DpOix, t)^ j d^c (c - u{x, tjff^ix, c, t) (4) 

Apparently, the above three relations correspond to 
the zero-th, first, and the trace of the second order 
hydrodynamic moments. It is well known that these 
three moments correspond to conservation laws and 
are invariant under any local collisions. 

Notice the density p is an overall multiplier on all 
moments, without loss of generality for the subse- 
quent analysis, we set it to unity. 

Now let us define an analogous hydrodynamic mo- 
ment expression in terms of summations over dis- 
crete velocity values below. 



M(") {x, t) = J2 Cac„ ■ • • Co /^^(a;, t) 



(5) 



2 



Or equivalently, in a Cartesian component form 

b 

M^,..,,^ {x, t)^Y. Co.,^,c^■^2 ■ ■ ■ Cc.,^Jl'{x, t) (6) 



In the above, we have assumed there are 6+1 num- 
ber of discrete _D-dimensional vector values in the 
basis discrete velocity set: {cq : a = 0, . . . , 6}. Sim- 
ilarly, we define an analogous equilibrium distribu- 
tion function. 



/^«((r,t) = iS„(0((r,t))exp 



where the macroscopic density, fluid velocity, and 
temperature are now defined in terms of moment 
summations instead. 



6 

q:=0 
6 



(8) 



Q=0 



In the above, Wa is a weighting factor that is at most 
dependent on 9{x^ t). Based on this fact, we can also 
re-express the discrete equilibrium distribution ([7]) 
in an alternative and simpler form: 



/^« = «)„(0)exp 
= WaiO) exp 



{Cg - Uf 

29 







exp 


261 



(9) 



by defining Wa{9) = Wa{9) exp[— Therefore, the 
discrete moment definition ^ can be re-expressed 
as, 



M. 



(n) 



r cg ■ ' 

V 6 



a=0 

X exp 







exp 


29 



(10) 



Having all the basic definitions above specified, we 
are now ready to prove several fundamental condi- 
tions for a lattice velocity basis supporting an n-th 
order hydrodynamic moment accuracy and its cor- 
responding form for the discrete equilibrium distri- 
bution function. These conditions are set forth for 



measuring any given lattice in terms of an intrinsic 
tensor: 



(n) 



(11) 



a=0 



Theorem 1 Discrete moment iW''-") is equal to the 
moment Af^") of the continuum Boltzmann kinetic 
theory, if the supporting lattice velocity basis satisfy 
the following conditions: 



e: 



(») 



gn/2A(n) 



0,2,4,. 



0, 



n = 1, 3, 5, . 



(12) 



In the above, A,-"'* ■ is the n-th order delta 
function defined as a summation of n/2 (n = 
even integer) products of simple Kronecker delta 
functions Si-^i.^ ■ ■ ■ and those from distinctive 

permutations of its sub- indices |ll|18|19|20j . There 
are (n — 1)!! (= (n — 1) • (n — 3) . . . 3 • 1) total num- 



ber of distinctive terms in A 

.(2) 



(n) 



For instance, 



A^;^ =(Sy,and 



ijkl 

ijklmn U kimn 



'■^ Imnj '■^ mnjh 



Sin A 



van J 
(4) 



(13) 



Obviously, a lattice velocity set that satisfy condi- 
tion (ini) for is 71-th order isotropic. 

Proof of Theorem 1: First we prove for the zero- 
th order moment, M.^'^^ — Af — 1. According to 
dH]) we have. 



: exp 



: exp 



u 

29 

29 



b 

E 

Q=0 



Wa exp 



1 



b 



(14) 



If (fT2|) is satisfied, then all odd valued I terms vanish, 
and the even valued terms become, 



= {2l-l)\Wu^ 



a=0 



21 



(15) 



In the above denotes a scalar product of two ten- 
sors. Therefore, (fH]) reduces to. 



3 



M(°) = cxp 



: exp 



20 
261 



1+E 



1^1 

oo 



1+E 



(2;-l)!!M2i' 

1 M^' 1 



(16) 



where the identity {21 ~ l)!!/(2/)! = 2^7'! is used. 
Since, 



1 V?' 



exp 



2^ 



(17) 



Substituting this into (|16p . we have proved that 
M(o) = 1. 

Next, we prove 

M(«) = m(") for n > 0. We 
start this by defining a partition function in discrete 
velocity space, 



Q = ^ Wq exp 



u 



exp 



2? 



(18) 



Notice the second equahty in the above is a result 
of the analysis of Af'^^ = 1. Consequently, we show 
that satisfying the second equality is a sufficient con- 
dition for achieving the correct hydrodynamic mo- 
ment for any integer n. First of all, we have the fol- 
lowing general relationship. 



It is easily shown that 
d^c(27r6i)-*e-«^ + ' 



exp 



'2Q 



Henceforth, we have shown that pT|) and ([20]) 
have exactly the same form. Subsequently, we have 
proved the theorem that M'^"^ — Af^"' for any 
positive integer n, if condition (|12p is satisfied. 

It is revealing to check a few obvious representa- 
tive examples. First of all, the first moment 



exp 



u 



du 



exp 



u 

20 



This is simply the fluid momentum or the fluid ve- 
locity. 

On the other hand, the second moment 



M(2 



exp 



u 

20 



exp 



20 



01 + uu 



where / is the second rank unity tensor. Hence the 
second moment has precisely the same form of the 
correct hydrodynamic momentum flux tensor. Fur- 
thermore, we have 



1 2 
-u 

2 



which is exactly the hydrodynamic total energy. 



■ exp 



: exp 



: exp 



u 

20 

20 
20 



Q=0 



ait" 



Wa exp 



a=0 



-Q 



(19) 



Since Q = exp [u'^/20] , then Eq. (fTO|) becomes 



M(") = exp 



u 

20 



du" 



(20) 



In comparison, from the continuum Boltzmann ki- 
netic theory, we have 



1 



dfc c • • • ccxp 



{2TT0)Dn 

•ri 

/ d^c c---c{2tt6 



{c~_uf 
20 



fin r 2 

-0n^ / i2n0)-^e-^ + ^ (21) 



3. Moment Accuracy for Lattices of Finite 
Isotropy 

In the previous section, we have proved that con- 
dition sufficiently ensures all moments defined 
via summations over discrete lattice velocity values 
are equal to that of the continuum Boltzmann ki- 
netic theory. However, such a condition is unneces- 
sarily too strong, because it requires the supporting 
lattice basis to have an infinite isotropy (i.e., n — s- 
oo). Obviously, no lattice velocity set containing a 
finite number of discrete values is able to meet such 
a requirement. Hence a realistic goal is to find a rela- 
tionship between the hydrodynamic moments up to 
a given finite order and the corresponding isotropy 
for the supporting lattice velocity basis. 

First of all, we notice the existence of a hierarchi- 
cal relationship among the hydrodynamic moments. 
Based on definition we have 



M(n) 



Hence, 



exp 



u 

20 



du 



Q 



(22) 



4 



M("+i) = exp 

,2 



: exp 



- exp 



w 

29 

26 
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2- 






20 




^ du) 
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du 
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if du 



9it 



exp 



20 



This gives the hierarchical relationship, 



(23) 



Using the hierarchical relationship ([23|) . all higher 
order moments are derivable starting from M^'^^ — 
1. More importantly, we realize that n-th order mo- 
ment Af^"' is an n-th order polynomial in terms of 
the power of the fluid velocity. That is, the highest 
power in M^^^ is it" . Since hydrodynamic moments 
up to a finite order only involve a finite power of fluid 
velocity, we expect moment accuracy up to a finite 
order can be achieved by a finite lattice set of ade- 
quate isotropy. Having established these properties, 
we arrive at the next theorem below. 

Theorem 2 // the supporting lattice velocity basis 
satisfies the following conditions 



E, 



(n) 



0, 



n = 0, 2, . . . , 2iV 

n — odd integer 



(24) 



and if the discrete equilibrium distribution function 
jei3,(Af) ^ truncation of the original exponential 
form by retaining terms only up to , then the dis- 
crete moment Af^") is accurate and equal to the mo- 
ment M^^^^ of the continuum Boltzmann kinetic the- 
ory for any n < N . N is any given finite positive 
integer. 

It is easily recognized that the basis lattice veloc- 
ity set satisfying the above condition must be 2N- 
order isotropic (c./., [Ill24j ). 

Proof of Theorem 2: We start by first examining 
the standard Maxwell-Boltzmann distribution ([3]), 
and express it in an expanded form in powers of fluid 
velocity u. This is very easily accomplished by tak- 
ing advantage of the following generating function 
for Hermite series, 



e.p[2tx-t^]=±^t 



(25) 



n=0 



where Hn[x) is the standard n-th order Hermite 
polynomial. Let us define the unity vector u = u/\u\ 



and |m| — yX^iii'"! magnitude, and ^ = 

c- u/ \/29. We can formally express the distribution 
® as. 



r{x,c,t)- 



(27r6 



(27r6 



- exp 


\ 








20 






exp 


c ■ u 







' 20 



' (:!,)« 

n=0 



{2tt0) 



A truncated series of the above can be de- 

fined by simply retaining the terms up to tt^. Based 
on the orthogonal property of the Hermite polyno- 
mials, namely 



dx e ^ Hm{x)Hn{x) — 0; Vm ^ n 



(27) 



and because c of power N can be fully represented 
by Hermite polynomials {Hn] n 0, . . . , N}, it is 
straightforward to see that moments up to A^-th or- 
der constructed out of /^'^ are identical to that of 
jeq,(N) ^ for the higher order terms in /'^'^ give van- 
ishing contributions due to orthogonality. 

Next, similar to the above, we expand the discrete 
distribution ([9|), and keeping terms only up to it^. 



re. 
J a 



U 



U 

20 



21 



Wa exp 

n=0 



n 



(28) 



where = Ca ■ it/ v20. Hence the task to prove 
Theorem 2 is to prove iVf'^"-' (Vn < N) generated 
by fa''-^^ is equal to M(") from the full Maxwell- 
Boltzmann distribution /^^^ or its truncation . 
According to definition ([T]) , we have 



M(") = J d^c c_^f'"'{x,c, t) 

n 

1 ^ 1 
^ I V- 



\ u 



X / d c c- ■ ■ c exp 



20 

£! 

20 



HniO (29) 



On the other hand, according to ©, we have 



5 



(JV) 



^From (|29|) and ((30|) , we see that both of these in- 
volve Hermite polynomials of orders no greater than 
N. Furthermore, a given Hermite function Hn{x) is 
a polynomial of a;™ (m = 0, . . . , < n). Therefore, 
both M(") and M(") involve powers of c (or Cq) 
from up to n + A^. Based this observation, we see 
that it is sufficient to prove M^"' = M(") (Vn < 
N) , if for all integer m < 2 the following property 
is satisfied. 



m 

Vto = 0, . . . , 27V 



= ^ Co • • • Cq (31) 



The result for the discrete summation is already 
given in the definition of and ([24|) . Hence it is 
suffice to just show that this is also true for the con- 
tinuum integration. In fact, according to the basic 
Gaussian integral property, we know that 



/im/2 A (m) 

0, 



d cexp 



21 



c 



Ci-, Ci. 



TO = 0,2,4, . 



,2N 



TO = 1,3,5, ...,2iV + l 



(32) 



Consequently, we have proved M^"^ = M^"^ (Vn < 
N), and thus Theorem 2. 

It is also worthwhile to note, without repeating 
the explicit steps of the above, that the same proof 
applies if the truncation of the exponential form 
is up to iV -f 1. Thus, we can retain an extra term in 
the expanded form. 

4. Discussion 

In this paper, we have presented and proved a 
set of fundamental conditions for formulating LBM 
models. Lattice velocity sets obeying these condi- 
tions automatically produce equilibrium moment 
accuracy to any given iV-th order. As demonstrated 
in [5], non-equilibrium moments are theoretically 
expressible as spatial and temporal derivatives of 
equilibrium moments. Therefore, achieving higher 



order moment accuracy enables accurate descrip- 
tion of fluid properties into deeper non-equilibrium 
regimes [21i22j . This is essential for physical prop- 
erties at finite Knudsen or Mach numbers that are 
beyond the Navier-Stokes representation. 

To make a more direct comparison with conven- 
tional LBM models, we rewrite in a more ex- 
plicit form (up to 0{u^)) below, 



Ca ■ U (Cq • m)2 



U 

29 



2612 

(Cq • u)^ {Ca ■ u)v? 

6^=^ 2^2 
(cq, • uY (cq • uf'u'^ 
24614 + ^ 

(Cq • m)^ [Ca-uY'V? {Ca ■ U)u'^ 



120615 



12614 



(33) 



It is immediately recognized that the series for 
most of the conventional LBM models terminate at 
0(^2) or 0{v?). For example, the so called D3Q15 
and D3Q19 correspond to the expansion up to 
0(^2) [15]. It can be directly verified that their 
underlying lattice velocity sets only satisfy the fun- 
damental conditions ((24| up to = 2, so that 
the higher order moment terms beyond 0(u^) can 
not be accurately supported. Furthermore, in these 
models, the temperature is fixed at = 1/3. An 
extended 34- velocity model exists |17l23j . and its 
temperature has a range of variation between 1/3 
to 2/3, and D3Q19 is its reduced limit as 6* = 1/3. 
But the moment accuracy is still N = 2. 

There are typically two approaches to construct 
lattice velocity sets obeying higher order of accura- 
cies {N > 2) according to (|24p. One approach is to 
rely on relations between discrete rotational symme- 
try and tensor isotropy |ll|24j . For instance, we can 
start with a lattice velocity set consisting of multi- 
ple lattice speeds, namely 

£ = £i U £2 • • • U Cm (34) 

where each of the subset is defined as 



^13 = {cq,/3; i = 0, 
/3 = 1,...,M 



All lattice velocities in each subset has the same 
magnitude, \ca,i3\ = cp. This way, the required 
isotropy can be imposed at each speed level. It 
has been shown that if such a velocity subset is 
parity invariant and obeys an n-th order isotropy 



6 



(n = even integer), then its basic moment tensor 
has the foUowing form [24 



E. 



_ u (-P-2)!! ^(„) 



(35) 



{D + n-2y.l 

and it vanishes for all the odd integer moments. Sub- 
sequently, we can assign a weighting factor fi;^(0) 
for each subset Cp,so that the overall condition (I24p 
is achieved by satisfying the following constraint on 
the weighting factors, 



M 

13=1 ^ 



{D-2)\ 



-n-2)\ 



(36) 



for n = 0, 2, ... , 2N. There are 27V + 1 such con- 
straints. Hence, it is necessary to include enough 
number of subsets and 'W/3{9) {(3 — I, . . . , M > N + 
1) in order to have a solution. Using such a proce- 
dure, a 59-velocity model in 3-dimension is formu- 
lated that satisfies p4|) up to iV = 3 with 6-th order 
tensor isotropy, so that the expansion in ([33]) can 
be carried to 0(u'*). Based on the analysis above 
and else where [5] , such an order of moment accu- 
racy is necessary for getting the correct energy flux 
in thermal hydrodynamics [25126127] . Another ap- 
proach is to form the discrete velocity sets via Gaus- 
sian quadrature for higher order models [5]. Indeed, 
P2p defines the precise requirement. The only dif- 
ference here is that the quadratures need to allow a 
variable temperature 0. This approach is relatively 
more straightforward, so that it enables a system- 
atic formulation of higher accurate LBM models to 
6-th, 8-th orders and beyond. There is also a simi- 
lar work recently by Sbragaglia et al on how to con- 
struct higher order isotropic moments [28] . 

The formulation described in this paper offers a 
rigorous measure for evaluating the order of accu- 
racy of a given LBM model. For future convenience, 
we may simply refer an LBM model that satisfies 
condition ^ to N-th order as "E(N)-accurate." 
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